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Abstract 

- - - In this paper, we introduce the notion of perfect space-time block codes (STBC). These codes have 

^^ ■ full rate, full diversity, non-vanishing constant minimum determinant for increasing spectral efficiency, 

^D i uniform average transmitted energy per antenna and good shaping. We present algebraic constructions 

, of perfect STBCs for 2, 3, 4 and 6 antennas. 
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Perfect Space Time Block Codes 



I. Introduction 

In order to achieve very high spectral efficiency over wireless channels, it is known that we need 
multiple antennas at both transmitter and receiver ends. We consider the coherent case where the receiver 
has perfect knowledge of all the channel coefficients. It has been shown [31] that the main code design 
criterion in this scenario is the rank criterion: the rank of the difference of two distinct codewords has 
to be maximal. If this property is satisfied, the codebook is said to be fully diverse. Once the difference 
has full rank, the product of its singular values is nonzero, and is defining the coding gain. Maximizing 
the coding gain is the second design criterion. Extensive work has been done on designing Space-Time 
codes that are fully diverse. 

We focus here on Linear Dispersion Space-Time Block Codes (LD-STBC), introduced in [21]. The 
idea of LD codes is to spread the information symbols over space and time. The linearity property of the 
LD-STBC enables the use of maximum likelihood (ML) sphere decoding [32], [20], which exploits the 
full performance of the code compared to other suboptimal decoders [8]. Consequently, research work 
has been done to construct LD-STBCs with more structure. One new property added has heen full rate, 
i.e., the number of transmitted signals corresponds to the number of information symbols to be sent, in 
order to maximize the throughput. In [9], it is shown how to construct full rate and fully diverse codes 
for the 2 transmit antennas case. This approach is generalized for any number M of transmit antennas in 
[11], [16]. A promising alternative approach based on division algebras is proposed in [28], where the 
authors construct non-full-rate and full-rate STBCs. A division algebra (as it will be detailed below) is 
an algebraic object that naturally yields a linear set of invertible matrices. It can thus be used to construct 
LD codes, since for any codeword the rank criterion is satisfied. 

In [4], [5], we have presented the Golden code, a 2 x 2 STBC obtained using a division algebra, which 
is full rate, full diversity, and has a nonzero lower bound on its coding gain, which does not depend on 
the constellation size. A code isomorphic to the Golden code was independently found by an analytical 
optimization in [33] and [10]. In [33, Theorem 1], it is also shown that, for 2 antennas, a sufficient 
condition for achieving the diversity-multiplexing gain frontier defined by Zheng and Tse [34] is exactly 
the lower bound on the coding gain. In [13], it has been shown in general that the nonzero lower bound 
on the coding gain is actually a sufficient condition to reach the frontier for any number of antennas. 

The goal of this work is to refine the code design criteria for LD-STBCs, asking for the three following 
properties: 

• A nonzero lower bound on the coding gain, which is independant of the spectral efficiency (non- 
vanishing determinant). 

• What we call a shaping constraint, to guarantee that the codes are energy efficient. 

• Uniform average transmitted energy per antenna is also required. 

We propose the so-called perfect codes that fulfill the above properties, and give explicit constructions 
in dimension 2, 3, 4 and 6 for 2 x 2, 3 x 3, 4 x 4 and 6x6 MIMO systems. 
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The paper is organized as follows. In Section |nl we detail the code design criteria and define precisely 
the notion of perfect codes. Since our code constructions are based on cyclic algebras, we begin Section 
irm by recalling how one can use cyclic division algebras to build fully-diverse and full-rate STBCs. 
We then explain further algebraic techniques useful to obtain the properties of the perfect codess. The 
following parts of the paper are dedicated to the code constructions. In Section Hvl we exhibit an infinite 
family of 2 x 2 perfect STBCs generalizing the Golden Code construction [5]. Then, we construct a 3 x 3, 
a4x4anda6x6 perfect STBC in Sections resp.lViUVl and 17111 

II. Problem statement 

We consider a coherent system over a flat fading M x N MIMO channel, where the receiver knows 
the channel coefficients (perfect CSI). The received matrix is 



NxT 



HAfxAf • XmxT + WatxT, (1) 



where X is the transmitted codeword of duration T taken from a STBC C, H is the channel matrix with 
i.i.d. Gaussian entries and W is the i.i.d. Gaussian noise matrix. Subscripts indicate the dimensions of 
the matrices. 

In this paper, we consider square (M = T) linear dispersion STBCs [21] with full-rate i.e., square 
codes with M degrees of freedom, using either QAM or HEX [14] information symbols. Since the 
codewords are square, we can reformulate the rank criterion saying that the codebook is fully diverse if 

I det(Xi - Xj)p / 0, Xi / Xj G C. 

By linearity, this simplifies to | det(X)p / 0, for all nonzero codeword X G C. 

Once a codebook is fully-diverse, the next step attempts to maximize the coding advantage, which is 
defined for LD-STBC by the minimum determinant of the code. We first consider infinite codes defined 
by assuming that the information symbols are allowed to take values in an infinite constellation. The 
minimum determinant of the infinite code Coo is 

<5min(Coo) = min |det(X)p. 

We denote by C the finite code obtained by restricting the information symbols to g-QAM constellations 
or gr-HEX. The minimum determinant of C is then 

<5„,in(C)= min |det(X)|2. 
o^xec 

In [28] as well as in all the previous constructions [11], [16], the emphasis is on having a non-zero 
minimum determinant. But since the minimum determinant is dependent on the spectral efficiency, it 
vanishes when the constellation size increases. 

Non-vanishing determinant. We say that a code has a non-vanishing determinant if, without power 
normalization, there is a lower bound on the minimum determinant that does not depend on the constel- 
lation size. In other words, we impose that the minimum determinant of the STBC is a constant Amin 
for a sufficiently high spectral efficiency. For low spectral efficiencies, it is lower-bounded by Amin. 
Non-vanishing determinants may be of interest, whenever we want to apply some outer block coded 
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modulation scheme, which usually entails a signal set expansion, if the spectral efficiency has to be 
preserved. 

A fixed minimum determinant is one of the two key properties of the perfect codes introduced in this 
work, the other one is related to the constellation shaping. 

Shaping, In order to optimize the energy efficiency of the codes, we introduce a shaping constraint 
on the signal constellation. It is enough to introduce this shaping constraint on each layer as the codes 
considered in this paper all use the layered structure of [12]. The g-QAM or g-HEX to be sent are 
normalized according to the power at the transmitter. However, since we use LD-STBCs, what is 
transmitted on each layer is not just information symbols but a linear combination of them, which 
may change the energy of the signal. Each layer can be written as Rv, where v is the vector containing 
the QAM or HEX information symbols, while R is a matrix that encodes the symbols into each layer. 
In order to get energy efficient codes, we ask the matrix R to be unitary. We will refer to this type of 
constellation shaping as cubic shaping, since a unitary matrix applied on a vector containing discrete 
values can be interpreted as generating points in a lattice. For example, if we use QAM symbols, we get 
the Z" (cubic) lattice. 

The last property of perfect codes is related to the energy per antenna. 

Uniform average energy transmitted per antenna. The ith antenna of the system will transmit 
the ith row of the codeword. We ask that on average, the norm of each row are similar, in order to 
have a balanced repartition of the energy at the transmitter. It was noticed in [28] that uniform average 
transmitted energy per antenna in all T time slots is required. 

We are now able to give the definition of a perfect STBC code. 

Definition 1: A square M x M STBC is called a perfect code if and only if: 

• It is a full rate linear dispersion code using M^ information symbols either QAM or HEX. 

• The minimum determinant of the infinite code is non zero (so that in particular the rank criterion is 
satisfied). 

• The energy required to send the linear combination of the information symbols on each layer is 
similar to the energy used for sending the symbols themselves (we do not increase the energy of 
the system in encoding the information symbols). 

• It induces uniform average transmitted energy per antenna in all T time slots, i.e., all the coded 
symbols in the code matrix have the same average energy. 

Let us illustrate the definition by showing that the Golden code, the 2x2 STBC presented in [5] is a 
perfect STBC. 

Example 1: A codeword X belonging to the Golden Code has the form 



\/5 



a{a + he) a{c + dO) 
ia{c + d6) a{a + h6) 



where a, 6, c, d are QAM symbols, = ^^^, 9 = ^-^, a = I + i - iO and a = l + i-i9. 

The code is full rate since it contains 4 information symbols, a, b, c, d. Let us now compute the minimum 
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determinant of the infinite code. Since aa = 2 + i, we have 



det(X) = ^^[(a + 60)(a + 6^)-i(c + d0)(c + d0)l 
5 



2-V 
By definition of a, 6, c, d, we have that the minimum of |a^ + ab — h"^ — i{cP + cd — (i^)P is 1, thus 



Jmm(Coo) = mill |det(X; 



,2_1 



OT^xec 5 

Thus the minimum determinant of the infinite code is bounded away from zero, as required. 
Let us now consider the diagonal layer of the code. It can be written 

1 la ad \ I a \ 
V5\ a ad J \ b J ' 

Since the matrix can be checked to be unitary, the cubic shaping is satisfied. 

Note in the second row of the codeword X the factor i, which guarantees uniform average transmitted 
energy since |i|^ = 1. 

This code has of course been designed to satisfy all the required properties. Its main structure comes 
from a division algebra, and the shaping is obtained by interpreting the signals on each layer as points 
in a lattice. In the following, we explain the algebraic tools we use, and show how to obtain codes with 
similar properties for a larger number of antennas. 

III. Cyclic algebras: a tool for space-time coding 

We start by recalling the most relevant concepts about cyclic algebras and how to use them to build 
full rate and fully diverse space-time block codes (see also [28] for more details ). We then explain how 
to add more structure on the algebra to get the other properties required to get perfect codes, namely, the 
shaping constraint and the non-vanishing determinant. We warn the reader that some algebraic background 
is required. If the reader is not familiar with the notions of norm, trace, Galois group, or discriminant, 
we recommand to read first the appendix U where these notions are recalled. 

A. Full rate and fully diverse STBCs 

In the following, we consider number field extensions K/F, where F denotes the base field. The set 
of non-zero elements of F (resp. K) is denoted by F* (resp. K*). 

Let K/F be a cyclic extension of degree n, with Galois group Ga.l{K/F) = (a), where a is the 
generator of the cyclic group. Let A = {K/F, a, 7) be its corresponding cyclic algebra of degree n, that 
is 

A=l-K®e-K®...® e^-^ ■ K 

with e e A such that le = ea{l) for all/ G K and e" = 7 G F*. 

Cyclic algebras provide families of matrices by associating to an element x ^ A the matrix of 
multiplication by x. 
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Example 2: For n = 2, we have A= 1- K @e- K with e^ = 7 and le = ea{l) for / G K. An element 
X ^ A can be written x = xq + exi. Let us compute the multiplication by x of any element y G A. 

xy = (xo + exi)(yo + eyi) 

= XQyQ + e(T{xQ)yi + exiyQ + -ia{xi)yi 

= koyo + 70-(2^i)yi] + e[cr(xo)yi + xiyo], 

since e^ = 7 and using the noncommutativity rule le = ea{l). 
In the basis {l,e}, this yields 

y xi o-(a:o) J \yi J 
There is thus a correspondance 

xo 70-(xi) \ 



In particular, 



X = Xo + exi (^ A ^^ 



eA 



xi a{xo) 




In the general case of degree n, we have for all Xk ^ K 



( Xk 

o-(xfe) 



Xk 



\ 



\ 







a 



ra-l 



and e 



{xk) I 



/ 
1 





7\ 




1 0/ 



Formally, one can associate a matrix to any element x ^ A using the map A^:, the multiplication by x 
of an element y ^ A: 

The matrix of the multiplication by A^^, with x = xq + exi + . . . + e"~^x„_i, is more generally given 
by 

/ Xo 7Cr(x„_i) 7Cr^(x„_2) ••• 7Cr"'"Ha^l) \ 
X\ cr(xo) 7Cr2(x„_i) ... 7Cr"~^(x2) 



7C7"-^(x„_i) 

^n— 1/ 



(2) 



x„_2 cr(a;„-3) o-^(x„„4) , 

y X„_i cr(x„_2) 0-2(x„_3) ... cr"~l(xo) j 

Thus, via A^;, we have a matrix representation of an element x ^ A. 

Let us show how encoding can be done. All the coefficients of such matrices are in K, K being a 
vector space of dimension n over F. Thus each Xj is a linear combination of n elements in F. The 
information symbols are thus chosen to be in F. If we consider QAM constellations with in-phase and 
quadrature ±1, ±3, . . ., the constellation can be seen as a subset of Z[i] := {a + 6z, a, 6 G Z} (Gaussian 
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integers). Since Z[z] C Q{i), we take F = Q(i) in order to transmit q-QAM. Similarly, in order to use 
HEX symbols, we see them as a subset of Z[j] := {a + bj, a,b ^'L} (Eisenstein integers) where j is a 
primitive 3rd root of unity {j^ = 1, j = e^*'^/^). We then take F = Q(j) with Z[j] C Q(j). Following 
the terminology of [28], we may say that the STBC Coo is over F. 
The following space-time block code is then obtained 

/ xq xi ... a;„_i \ 

7cr(x„_i) a{xo) ... a{xn-2) 



Cr, 



Xj £ K, i = 0, 



,n 



(3) 



[\ja^~\xi) 7a"-i(x2) ... ct"-1(xo) / 

Since each codeword X contains n coefficients Xj, each of them being a linear combination of n 
information symbols, cyclic algebras naturally yields full rate LD-STBCs. 

Definition 2: The determinant of the matrix Q (which is also the determinant of a codeword Q) is 
called the reduced norm of x, a; G A. 

The key point of this algebraic scheme is that we have a criterion to decide whether the STBC Coo 
satisfies the rank criterion. Namely, when the cyclic algebra is a division algebra, all its elements are 
invertible; hence the codeword matrices have non zero determinants. 

Proposition 1: [28] The algebra A = {K/F, a, 7) of degree n is a division algebra if the smallest 
positive integer t such that 7* is the norm of some element in K* is n. 

So at this point, by choosing an element 7 such that its powers are not a norm, the codebook Coo defined 
in (|3ll is a fully diverse LD-STBC with full rate. 



B. The shaping constraint using complex algebraic lattices 

The shaping constraint requires that each layer of the codeword is of the form Rv, where R is a 
unitary matrix and v is a vector containing the information symbols. Let K = F{6) and {1,9, ... , 9"'^^} 
be a F-basis of K. Each layer of a codeword X as in (|3ll is of the form 



/l 
1 



1 a 



a{9) 
""i(^) 



nn— 1 



_7i— 1 1 Qri—l\ 



\ 



J 



( 



Ul,0 



\ 



( 



\ Ul-a-l / 



XI 



\ 



(4) 



V ^"-'(^0 / 



for xi = X]fc=o '^i,k(^^- Since ui^^ takes discrete values, we can see the above matrix multiplication as 
generating points in a lattice. The matrix R is thus the generator matrix of the lattice, and the lattice 
obtained is given by RR^, its Gram matrix. We would like R to be unitary, which translates into saying 
that the lattice we would like to obtain for each layer is a Z[i] "-lattice, resp. a Z[j]"^-lattice, since QAM 
and HEX symbols as finite subsets of Z[z], resp. Z[j]. Note that the matrix R may be viewed as a 
precoding matrix applied to the information symbols. 

Finally, note that the 2n^-dimensional real lattice generated by the vectorized codewords where real 



and imaginary components are separated, is either 



72n^ 



(for QAM constellation) or A^ (for HEX 
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constellation), where A2 is the hexagonal lattice [7], with generator matrix 




Interpreting the unitary matrix R as the generator matrix of a lattice allows us to use the well studied 
theory of algebraic lattices [1], [2], [24]. The key idea is that the matrix R given in (HJ needs to contain 
the embeddings of a basis, but this basis does not need to be a basis of the field K. It can be a basis of 
a subset of K, and in fact it will be a basis of an ideal of K. 

Let K he. Si Galois extension of F = Q(i) (resp. F = Q(j)) of degree n, and denote by Ok its ring 
of integers. Let Q(0) be a totally real Galois number field of degree n with discriminant coprime to the 
one of F, that is {dF,dQ/g\) = 1. In the following, we focus on the case where K is the compositum 
of F and Q(^) (that is, the smallest field that contains both). We write the compositum as i^ = FQ{6) 
(see Fig. [l}- This assumption has the convenient consequence that [30, p. 48] 

d-K = dn(o)dp, (5) 



where dp = -4 for F = Q(i) and dp = -3 for F -- 

Denote by {0"^}^^^ the Galois group Gal(i^/F). 

Definition 3: We denote by A^(T) the complex algebraic lattice corresponding to an ideal T C Ok 
obtained by the complex embedding a of K into C" defined as 

cr : i^ ^ C" 

X ^ a{x) = (<Ti(x), . . . , an{x)) 
The basis of A'^(X) is obtained by embedding the basis {z^A,}fc=iOf ^- Consequently its generator matrix 
is similar to the matrix R in Q, where the basis {1,0, ... , ^"~^} is replaced by the ideal basis {i'k}^=i- 
Its Gram matrix G is thus given by 

G = {TrK/F{^k'A))k,i=i, 

where x denotes the complex conjugation of x. When F = Q(j), since Gal(Q(j)/Q)=(r), with t{j) = 
j'^ = j, we have that r coincides with the complex conjugation. 

We explain now how to choose an ideal T C Ok in order to get the rotated versions of the Z^" or A2 
lattices. First consider the real lattice A(T) obtained from A'^(T) by vectorizing the real and imaginary 
parts of the complex lattice vectors. We want A(T) to be a rotated version of Z^" or ^2 • The basic idea 
is that the norm of the ideal X is closely related to the volume of A(X). We will thus look for an ideal 
with the "right" norm. 

• Consider the ramification in K/Q, that is the way prime numbers in Z may factorize when considered 
in K (for example, 5 is prime in Z but is not prime anymore in Q(i), since 5 = (2 + i)(2 — i)). We 
say that a prime p^ ramifies if {pk)OK = Yie-^M where e^ > 1 [25, p. 86] for some k (or in words, 
the primes which when factorized in K have factors with a power greater or equal to 2). The prime 
factorization of the discriminant dK/q = Y[Pk' contains the primes pk which ramify [25, p. 88]. 
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Considering real algebraic lattices K{Ok) [1], we know that vo\{K{Ok)) = '^^"' \ \dK/Q\- We 
look for a sublattice A(X) of K{Ok), which could be a scaled version of Z^" (resp. A2), i.e., 
A(T) = (y^Z)^" (resp. {cA2)'^) for some integer c. 
Since A(J) is a sublattice of A(C'i<-), vol(A(C'i^)) = 2""Y^K?d must divide 

'' VOl((^Z)2") 

vol((c.42)") = C" 






vol(A(X)) 



i.e., d^/Q = n^fc" divides 22"c2" (resp. S^c^"). 

This gives a necessary condition for the choice of T. In terms of norm of the ideal T [25, p. 69], 

we need 



vol(A(J)) 



vo1(A(Ok)) 



(2c)" 



(V3c)" 



N{1) = \Ok/T\ 

Recall from ^ that cLk = 2^" • d'L,g., when K is the compositum of 
discriminants and that d/^ = 3" • d'L,^^, when K is the compositum of 
discriminants. 



and 
and 



(6) 

with coprime 
with coprime 

n-rfc/2. 



• In order to satisfy Q, we must find an ideal X with norm np,=^2^fc ^'' (resp. rips^aPfc ^'^ )■ 
This procedure helps us in guessing what is the "right" ideal T to take in order to build a I?^ or A2 
lattice. To prove that we indeed found the "right" lattice it is sufficient to show that 

^"^K/pinvj) = 6ij, i,j = l,...,n (7) 

where {i^i}f^i denotes the basis of the ideal X, and 6ij is the Kronecker delta. 

Note that the lattice does not exist on all field extensions K/F. Once we have a cyclic field extension 
where the lattice exists, we define a fully diverse full rate codebook which furthermore satisfies the 
shaping constraint as 

• • • Cr{Xn~2) 



Cx 



' ( Xq 

7Cr(x„_i) 



cr(xo) 



[ V 7a"-i(xi) ^a^~^{x2) 



a 



n~l 



a^i G X C Ok, i = 0, . . . ,n — 1 



(^0) / 



> . 



(8) 



C. Discreteness of the determinants 

The goal of this section is to show how to get codes built over a cyclic algebra A = {K/F, a, 7) so 
that their determinants are discrete. One condition will appear to be 7 G Op, the ring of integer of F. 
This contrasts with the approach of Sethuraman et al. [28, Proposition 12], where the element 7 was 
chosen to be transcendental. Hence, the cyclic division algebra {K{'^)/F{'^),a, 7) is used, which ensures 
that the minimum determinant is non-zero. Unfortunately, this approach yields a vanishing minimum 
determinant, when the constellation size increases. 

In [3], it has been shown for 2 x 2 STBCs, by an explicit computation of the determinant, that the 
reduced norm of the algebra (see Def. ^ is linked to the algebraic norm of elements in K. Since the 
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norm of an element in K belongs to F, restricting the codeword matrix elements to be in Ok and taking 
7 G Op then gives discrete values of the determinants for the codewords of 2 x 2 STBCs. The same 
result has also been used in [5], for the Golden code. However, an explicit determinant computation is 
no more possible in higher dimensions. We thus invoke a general result that guarantees the reduced norm 
to be in F. 

Theorem 1: [26, p. 296 and p. 316] Let A = (i^/F, cr, 7) be a cyclic algebra, then its reduced norm 
belongs to F. 

Corollary 1: Let A = {K/F, a, 7) be a cyclic algebra with 7 G Op. Denote its basis by {1, e, . . . , e"~^}. 
Let X G ^ be of the form 

X = xq + exi + . . . + e^~^Xn-i 

where Xk G Ok, A; = 0, . . . , n — 1. Then, the reduced norm of x belongs to Op- 
Proof: Recall from Definition |2l that the reduced norm of x is the determinant of its matrix 
representation. Since Xk G Ok implies a{xk) G Ok for all k and 7 G Op by hypothesis, all coefficients 
of the matrix representation belong to Ok, hence so does its determinant. By Theorem^ the reduced 
norm of x belongs to F, so it belongs to Ok n F = Op. ■ 

Corollary 2: The minimum determinant of the infinite code with T = Ok defined in ^ is 

^mm(Cc)^) = 1. 

Proof: Since we only consider Op = TL\}\ (resp. Op = Z[j]), the determinants of the codewords 
form a discrete subset of C: 

det(X) G Z[i] (resp. G Z[j]). 

Then 5min(Cc)^.) = minx^o I det(X)p = 1 as the minimum is achieved by taking the codeword with 
xo = 1 and x^ = for A; = 1 . . . n — 1, corresponding to a single information symbol uqq = 1 and all 
the remaining v? — 1 equal to 0. ■ 

Let us give as example the 3x3 case to show that things become more complicated than the 2x2 case 
when the dimension increases, so that the general Theorem^ is required. 

Example 3: Consider a cyclic algebra A = {K/F, a, 7) of degree 3 with 7 G Op. Let x = xq + exi + 
6^X2, which can be represented as 



X 



Xo Xi X2 

7a(x2) cr(xo) cr(xi) 
7ct2(xi) 7C72(x2) cr2(xo) 



The norm of x is given by the determinant of X: 

det(X) =7^X2Ct(x2)ct2(x2)+ 

7{-XoCr(xi)o-^(x2) - 0-(xo)cr^(xi)x2 - cr^(xo)xicr(x2) + Xicr(xi)o-^(xi)} 
+ XoCr(xo)cr^(xo) 
=iV(xo) + -fN{xi) + 7^iV(x2) - 7Tr[xocj(xi)cj2(x2)] 
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Obviously the norm of the algebra is not only related to the norm of the number field, as in dimension 
2 where det(X) = N{xi) - 'yN{x2) [5]. 

When considering a particular case, it is possible to conclude that det(X) still belongs to F, either as 
in Example |3] by finding an expression in terms of norms and traces, or by noticing that the determinant 
is invariant under the action of a. Since the expression in larger dimensions gets more complicated, for 
the general case, we simply use Theorem^ 

Note that at this point, we have all the ingredients to build perfect codes. Assume there exists 7 G Op 
such that none of its powers is a norm. Then the code Cj defined in ^ is fully diverse and full rate, it 
has the required shaping constraint, and we have just shown that its determinant is discrete. In order to 
conclude, it is now enough to take I7P = 1, to guarantee uniform average transmitted energy per antennas. 
Before summarizing our approach, we now give an explicit bound on the minimum determinant. 



D. The minimum determinant 

We discuss now the value of the minimum determinant of the codes. Depending on whether the ideal 
T introduced in subsection IIII-EI is principal (i.e., generated by one element), we distinguish two cases. 
We show that if T is principal, then the minimum determinant of the infinite space-time code Cj is easily 
computed. Otherwise, we give a lower bound on (^min(Cj). 

Let us first assume T = {a)OK is a principal ideal of Ok- For all x G T, we have x = ay for some 
y G Ok- Notice that in this case, codewords are of the form 



X 



a 
a{a) 

••• 











a 



n-l 



(«) 



yo 
i(y{yn-i) 



yi 
cr(yo) 



yn-i 
o-(yn-2) 



7cr 



n-l 



(yi 



7cr 



n-l 



(y2 



a 



n — 1 



(yo 



(9) 



where yi G Ok, i = 0,...,n — 1. Since 7 G Op, the determinant of the second matrix is in Op and by 
Corollary |2l its square modulus is at least 1. We deduce, recalling that F = Q(i) or Q(j), that 



^(^^) = ni^i^. |det(X)|2 = \NK/F{a)\'' = NK/qia). 



(10) 



The last equality is true since the complex conjugation is the Galois group of F/Q. Thus 

n— 1 n— 1 

fc=0 fc=0 

where a^, and a^, k = 0, . . . ,n — 1, give the 2n elements of the Galois group of K/Q. 

Since K is the compositum of F and a totally real field Q(0) and we require the cubic shaping, we 
can go a little further. 

Proposition 2: Let Cj be a perfect code built over the cyclic division algebra A = {K/F, a, 7) of 
degree n where 7 G Op, K = FQ{6) and I is principal. Then 

1 



■'minv'-'X 



(Ci 



d, 
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where (iQ(6)) is the absolute discriminant of 

Proof: Let {iyi}f^i be a basis of the principal ideal I = {a)OK and A(X) denote the real lattice 
over Z. Recall [1] that 

det(A(J)) = vol(A(T))2 = A-'^NilfdK (11) 

where dx denotes the absolute discriminant of K. Using ^ and considering the real lattice, we have 

for F = Q(i) 

det(z2") = 1 = 4-"Ar^/Q(a)24(,)4", 

and for F = Q(j) 

det(A^) = (3/4)" = 4-"Ar^/Q(a)24(,)3". 

Both cases reduce to 

(a) 



and we conclude using (ITOt . ■ 

We consider now the more general case, where we make no assumption on whether T is principal. We 

have the following result. 

Proposition 3: Let Cj be a perfect code built over the cyclic division algebra A = {K/F, a, 7) of 

degree n where 7 G Op- Then 

6min{Cl) G N{T)Z, 

where iV(X) denotes the norm of T. 
Proof: Recall first that 

n 

det(X)= J];sgn(s)J](X)fc,,(fc), 

sG5„ k=l 

where Sn is the group of permutations of n elements, and sgn denotes the signature of the permutation. 
Denote by Z"^ the action of the Galois group on Z. Since (X)fc^(^,) G Z"^ for all k, we get [15, p. 118] 

det (X) G n ^''=^K/F {!) Ok , 

cr&Ga\{K/F) 

where Mx/pi^) stands for an ideal of Op called the relative norm of the ideal Z. The notation NiZ) 
emphasizes the fact that in the case of the relative norm of an ideal, we deal with an ideal, and not with 
a scalar, as it is the case for the absolute norm N{Z) of an ideal. 

Note that the above equation means that det(X) belongs to an ideal of Ok- By Corollary [Q we deduce 
that 

det(X) G Opr\MK/p{l)OK=MK/p{l), 

which means that det(X) is actually in an ideal of Op. Thus | det(X)p G J\fp/iQ{J\fK/p{Z)), since again 
F = Q(i) or Q{j). We conclude using the transitivity of the norm [15, p. 99] 

min |det(X)|2 G Mk/q{Z) = N(Z)Z. 

-A.GCx;-A.7=U 
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Bounds on (5min(Cj) are easily derived from the above proposition. 

Corollary 3: Let Cj be a perfect code built over the cyclic division algebra A = {K/F, a, 7) of degree 
n where 'j e Op and K = FQ{e). Then 

N{I) = < (5min(Cj) < — — — — — rninA^;^/Q(x) 

dQ(e) vol(Ac(J)) xei '^ 

Proof: The lower bound is immediate from Proposition|3]and the equality comes from dlU . similarly 
as in the proof of Proposition |2l 

An upper bound can be obtained as follows. We take xq 7^ G X, xi = . . . = x„„i = 0, which 
yields as determinant Nj^/q^xq). Thus mindet(X) = miux^i Nj^/q{x). Since the ideal Z may give a 
scaled version of the lattice Z[i]" (resp. Z[j]"), a normalizing factor given by the volume of the lattice 
is necessary to make sure we consider a lattice with volume 1. ■ 

The result obtained in (fTOb for the principal case alternatively follows: 

Corollary 4: If T = {a)OK is principal, then 

Proof: If T is principal, the lower and upper bounds in Corollary |3l coincide. ■ 

E. Summary of our approach 

Let us summarize the techniques explained above, and give the steps we will follow in the next sections 
to construct perfect codes: 

1) We consider QAM or HEX symbols with arbitrary spectral efficiency. Since these constellations 
can be seen as finite subsets of the ring of integers Op = Z[i] (resp. Op = Z[i]), we take as base 
field F = Q(i) (resp. F = Q(j)). 

2) We take a cyclic extension K/F of degree n = M with Galois group Gal{K/F) = {a) and build 
the corresponding cycUc algebra: 

A = {K/F,a,^). 

We choose 7 such that I7I = 1 in order to satisfy the constraint on the uniform average transmitted 
energy per antenna. 

3) In order to obtain non-vanishing determinants, we choose 7 in Z[i], resp. in Z[j] (see Sec. Illl-Ct . 
Adding the previous constraint I7I = 1, we are limited to 7 G {l,i,— 1,— i} C Z[i] or 7 € 
{1, J, j^ -1, -J, -f} C Z[j], respectively. 

4) Among all elements of A, we consider the discrete set of codewords of the form x = xq + xie + 
. . . + x„_ie"~^, where Xj G Z, an ideal of Ok, the ring of integers of K. This restriction on the 
coefficients guarantees a discrete minimum determinant (see Section Illl-Ct . We thus get a STBC 
of the form 

/ Xo Xi ... X„_i \ 

7Cr(x„_i) (t(xo) ... cr(x„_2) 



Cx 



a^i G 2r C Ok, i = 0, ■ ■ ■ ,n — 1 



[V7CT"-i(xi) ja^-Hx2) ... a"-i(xo) / 



(12) 
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The n? information symbols U£ ^ G Op are encoded into codewords by 

n-l 

^1 = ^^ Ui,k'^k i = 0,...,n-l 
fc=o 

where {i'k}^Zo ^^ ^ basis of the ideal I. 

5) We make sure to choose an ideal I C Ok so that the signal constellation on each layer is a finite 
subset of the rotated versions of the lattices Z^" or Ag . 

6) We show that A = {K/F, a, 7) is a division algebra by selecting the right 7 among the possible 
choices which reduces to show that 7, . . . , 7""^ are not a norm in K*. 

Since the desired lattice does not always exist, we need to choose an appropriate field extension K that 
gives both the lattice and a division algebra. Note that, in building a cyclic algebra A = {K/F, a, 7) for 
STBCs, the choice of 7 is critical since it determines whether ^ is a division algebra. It is furthermore 
constrained by the requirement that I7I = 1, so that the average transmitted energy by each antenna in 
all time slots is equalized, and to be in Op to ensure the discreteness of the determinant. 

Remark 1: The construction of the codes involves a lot of computations in number fields. Some of 
them are done by hand, some of them are computed with the computational algebraic software Kant [35]. 

IV. An infinite family of codes for 2x2 MIMO 

In this section, we generalize the construction given in [5] to an infinite family of codes for 2 x 2 
MIMO. 

Let p be a prime. Let K/Q{i) be a relative extension of degree 2 of Q(i) of the form K = Q(i, ^/p). 
We can represent K as a vector space over Q{i): 

K = {a + hy/p I a,&eQ(i)}. 

Its Galois group Gal(K/Q(z)) = {a) is generated by a : ^Jp h^ — ^/P- The corresponding cyclic algebra 
of degree 2 is ^ = (/s:/Q(i), a, 7). 

We prove here that when p = 5 (mod 8), 'j = i, and using a suitable ideal I C Ok, we obtain perfect 
codes following the scheme of Sec. IIII-EI 

A. The lattice Z[i]'^ 

We first search for an ideal I C Ok giving the rotated Z[i]^ lattice. We use the fact that Z[i]^ is the 
only unimodular Z[i]-lattice in dimension 2 [27]. Hence it is enough to find an ideal X such that A'^(X) 
is unimodular. By definition, a unimodular lattice coincides with its dual defined as follows. Let A'^(T) 
be a complex algebraic lattice with basis {vi, V2} = {a{vi),cr{v2)} following the notations of Section 

mrB] 

Definition 4: The dual lattice of A^(X) is defined by 

A^(J)# = {x = aivi + a2V2, ai, 03 G Q(i)| (x, y) G Z[i], Vy G A'=(J)} 
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where the scalar product between the two vectors can be related to the trace of the coiTcsponding algebraic 
numbers as 

(x, y) = Tr/^/Q(i) (xy) . 
The dual of a complex algebraic lattice can be computed expUcitly. Recall that the codifferent [30, p. 
44],[1] is defined as 

D-i^ = {xe K \\/ae OK,Trii/F{xa) G Op}. 

Lemma 1: We have A^I)* = A= (J#) with 



where T^x]n(i) denotes the codifferent (defined above). 

Proof: Let x G l~'^'D~^\Q,^y For all y e 1, v^e have to show that TTCK/Q{i){xy) ^ ^W- Since 
X = uv, with u G I^^ and v G '^x)n(i)' we have xy = uyv, with uy G Ok- The result follows now 
from the definition of VjJ'ir^r-y ■ 



Let K = Q{i, ^/p), with p = \ (mod 4). The factorization of p in Ok is [6] 



(p)Oi^=j2.T' (13) 

where T, Z are prime conjugate ideals. 

Proposition 4: The Z[i]-lattice -^A'^(T) is unimodular. 

Proo/.- Note first that Vkj^kC) = ^Q(v^)/Q = i^/P)^<Q>{^) = (Vp)- Using Lemma [l] and O, we 
compute the dual of X, 

1 



Now the dual lattice is 



p 

-^A^(J)) = v^ (a^(J)#) = -^A^(J) 
which concludes the proof. ■ 

B. The norm condition 

The last step is to prove that the algebra A = {K/Q{i),a, i) is a division algebra. In order to do that, 
we have to show (see Proposition Qi that 7 = i is not a norm in K/Q{i). 

We first recall the characterization of a square in finite fields. Let p be a prime and denote by GF{p) 
the finite field with p elements. 

Proposition 5: Let x G GF{p)*. We have 

X is a square <^=^ x^~ = 1. 
Proof: See [22]. ■ 

Corollary 5: If p = 1 (mod 4) , — 1 is a square in GF{p). 
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Let US come back to our case where p is a prime such that p = b (mod 8) and K = Q(i, ^/p) is a relative 
extension of Q(i). Let x = a + hy/p G K, a,b e Q{i). Its relative norm is 

NK/Q{i){x) = (a + b^){a - b^) = a^ - pb^ . (14) 

Our goal is to show that the equation Ni^iqij\{x) = i has no solution. As in [5], we prove that this 
equation has no solution in the field of p-adic numbers Qp, and thus, no solution for x ^ K. Let 
Zp = {x G Qp|fp(a;) > 0} be the valuation ring of Qp, where Vp{x) denotes the valuation of a; in p (that 
is, the power at which p appears in the factorization of x). First, we check that i G Zp. In fact, there are 
embeddings of Q(i) into Qp if X'^ + 1, the minimal polynomial of i, has roots in Zp. Using Hensel's 
Lemma [17, p. 75], it is enough to check that -1 is a square in GF{p). By assumption, p = 5 (mod 8), 
thus p = I (mod 4), then, by Corollary 13 -1 is a square in GF{p). 

Proposition 6: The unit i G Z[i] is not a relative norm, i.e., there is no x G K such that Ni^/Q(^i^{x) = i 
where K = Q(a/P, i) with p = 5 (mod 8). 

Proof: This is equivalent, by (fT4l) . to prove that 

a^ -pb^ = i, a,6G Q(i) (15) 

has no solution. Using the embedding of Q(i) into Qp, this equation can be seen in Qp as follows: 

a^ -pb^ = y + px, a,b e Qp, x,y e Zp, (16) 

where y^ = —1. If there is a solution to (I15t . then this solution still holds in Qp. Thus proving that no 
solution of (fT6b exists would conclude the proof. We first show that in (I16t . a and b are in fact in Zp. 
In terms of valuation, we have 

i/p(a^ - pb"^) = yp{y + px). 

Since x G Zp and y is a unit, the right term yields Up{y + px) > inf{fp(y), i/p(x) + 1} = 0, and 
we have equality since the valuations are distinct. Now the left term becomes = Up{a'^ — pb"^) = 
inf{2fp(a), 2vp{b) + 1}. The only possible case is Vp{a) = 0, implying a ^"Lp and consequently b G Zp. 
We conclude showing that 

a — pb =y+px, a,b,x,y (^ Zp (17) 

has no solution. Reducing (mod pZp), we see that y has to be a square in GF{p). Since y^ = —1, 
y(p-i)/2 = (|_i)(p-i)/4 _ _-Y ]^y choice of p = 5 (mod 8). By Proposition |5j y is not a square, which 
is a contradiction. ■ 

Remark 1: This result does not hold for p = 1 (mod 8) since, in this case, y(P~i)/2 = (_i)(p-i)/4 = \ 
and we get no contradiction. The fact that this proof does not work anymore is not enough to restrict 
ourselves to the case p = 5 (mod 8). We thus give a counterexample. 

Example 4: Consider K = Q (\/l7, i) , and x = ^^ ' — ^~ j — . It is easy to check that Nj^/q^^^^ (x) = 
i. 
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C. The minimum determinant 

We first show that the ideal T in (I13t is principal for all p = 1 (mod 4). Since N{T) = p, it is enough 
to show that there exists an element a G X with absolute norm N^/Qia) = p. Using the fact that 
p = u^ + v"^ for some n, u £ Z (that can be computed) [25], the element a = \/n + iv has the right 
norm and generates I (resp. a = \/u — iv generates T). Now, take 9 = — ^ and let 6 = — ^ be its 
conjugate. We have Ok = "^[9]. The codewords have the form 






a{a + h9) a{c + d9) 
ia{c + d9) a{a + b9) 



with a,b,c,d € Z[z]. Each layer of the STBC can be encoded by multiplying the vectors (a,b)'^ and 
(c, d)^ by the matrix 

a a9 

a a9 

which generates the Z[i]^ lattice. We observe that this lattice generator matrix may require basis reduction 
in order to be unitary. 

Determinants are given by 

det(X) = -NK/Q(i)ia) {NK/Q(i){a + b9) - iNK/Q(^i)[c + d9)) . (18) 

As the second term in ( fTHl only takes values in Z[i] and its minimum modulus is equal to 1 (take for 
example o = 1 and b = c = d = 0), we conclude that 

SminiCoc) = — |^E-/Q(j)(a)P = —NK/qia) = - ■ (19) 

Remark 2: As p = 5 (mod 8), the largest minimum determinant is given by p = 5 corresponding to 
the Golden code [5]. 

V. 4 X 4 PERFECT STBC CONSTRUCTION 

As for the 2x2 case, we consider the transmission of QAM symbols, thus, the base field is F = Q(i). 

Let 6* = Ci5 + Cili^ = 2cos(ff ) and K be Q(i, 9), the compositum of Q{i) and Q{9). Since (p{15) = 8 

{ip is the Euler Totient function), [Q{9) : Q] = 4, and thus [Q(i, 9) : Q(i)] = 4. The discriminant of 

is dQfg\ = 1125 and the minimal polynomial P0{X) = X^ — X^ — AX"^ + 4X + 1. The extension 

(i) is cyclic with generator cr : Ci5 + Cis"*^ ^^ C15 + Ci5^- 

The corresponding cyclic algebra of degree 4 is ^ = (K/Q(i),cr, 7), that is 

A = l-K®e-K®e'^-Kee^-K 

with e G K such that e^ = 7 G F* and le = ea{l) for all / G i^. In order to obtain a perfect code, we 
choose 7 = i. 
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A. The 'L[iY' complex lattice 

We search for a complex rotated lattice Z[i]^ following the approach given in liii-BI Since the relative 
discriminant of K is dj^ip = dqrg^ = 1125 = 3^ • 5^, a necessary condition to obtain a rotated version 
of Z[z]^ is that there exists an ideal T C Ok with norm 45 = 3^ • 5. The geometrical intuition is that the 
sublattice A(X) has fundamental volume equals to 2~^^/dKN{Z) = 3^ • 5^ = VTs , which suggests that 
the fundamental parallelotope of the algebraic lattice A(X) could be a hypercube of edge length equal to 
^/T5. 

An ideal T of norm 45 can be found from the following ideal factorizations 



(3)Ox 



^3^3 



{h)OK = XlXs' 

Let us consider X = X3 • X5. It is a principal ideal X = (a) generated by a = (1 — J>i) + iO'^. 

A Z[i]-basis of (a) is given by {a0*}?^Q. Using the change of basis given by the following matrix 

/ 1 \ 



V 






1 





-3 1 


-1 


-3 1 1 



/ 



one gets a new Z[z]-basis 

{vk]l=i = {(1 - 3i) + io'^, (1 - ^i)e + ie^ -i + (-3 + 4i)0 + (1 

Then by straightforward computation we can check that 

(i^ki^i) = Ski k,£ = l, 
9, Tr, 



{-i + i)-39 + e^ + e^}. 



Y^Trx/Q(i)i 



using TrQ(e)/Q(6') 

compute the diagonal coefficients, 



1> TrQ(e)/Q(6'^) = 9, TrQ(5i)/Q(6'^) = 1, TrQ(e)/Q(6''') 



29. For example, we 



Tr/f/Q(i)(|i^fc| ) 



TrK/Q(i){W - 66^ + 9'^) = 15 if A; = 1 

Tri^/Q{i) (1 + 39 + 9^-6^) = 15 if k = 2 

TrK/Qii)i5 + 69-9^ 
Tri^/Q»(-5^ + 2e2 



-29^) 
29^) -- 



: 15 if A: 
15 if k 



The unitary generator matrix of the lattice is given by 

1 



V 15 

/ 0.2582 - 0.3122i 
0.2582 + 0.0873i 
0.2582 + 0.2136i 

y 0.2582 - 0.7633i 



0.3455 - 0.4178i 

0.4718 + 0.1596i 

-0.5051 - 0.4178i 

-0.054 + 0.1596i 



-0.4178 + 0.505H 
0.1596 + 0.054i 

-0.4178 - 0.3455i 
0.1596 - 0.4718i 



-0.2136 + 0.2582i \ 
0.7633 + 0.2582i 
0.3122 + 0.2582i 

-0.0873 + 0.2582i y 
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B. The norm condition 

We now show that A = {K/Q{i),a,i) is a division algebra. By Proposition Q we have to check that 
ibi and —1 are not norms of elements in K. 
Lemma 2: We have the following field extensions: 

q{i)cQ{i,V5)cK 

Proof: We show that Q(z, \/5) is the subfield fixed by (o-^), the subgroup of order 2 oiGa[{K/Q{i)) = 
(a). Let a^ : Ci5 + Cil,^ ^ Ci^s + Cfs^ and x = Y.l=oak{C + C"^)^ Ofc e Q{i), be an element 
of K. It is a straightforward computation to show that a'^{x) = x implies that x is of the form 

x = ao + a3(Ci5 + Cii) = ao + ^3^^^ G Q(«, V5). ■ 

Proposition 7: The algebra A = {K/Q{i),a,i) is a division algebra. 

Proof: We start by proving by contradiction that iti are not a norm. Suppose ±i is a norm in K*, 
i.e., there exists x G ii'* such that Nj^/Qi^i-^{x) = ±z. By Lemma Sand transitivity of the norm, we have 

NKmi)ix) = A^Q(i,V5)/Qw(^i^/Q{i,V5)(^)) = ±^- 

Thus ±i has to be a norm in Q{i, \/5)- By Proposition 15] in the case p = 5, we know i is not a norm. In 
order to show that —i is not a norm, it is enough to slightly modify the proof of Proposition |6l Eq. (fTTl 
becomes, with p = 5, 

a — 5b = —y + 5x, a,b,x,y (^ Z^. 

Reducing (mod 5), we see that in order for this equation to have a solution, y has to be square in 
GF{5). Since {—y)^^^'^''^ = (— y)^ = 'iP' = —1. y cannot be a square (see Proposition [Sjl and we get a 
contradiction. 

The previous argument does not apply for —1 since it is clearly a norm in Q(i, \/5)/Q(i). The proof 
that —1 is not a norm uses techniques from Class field theory and is given in Appendix II VI ■ 

C. The minimum determinant 

From dlOl i. or similarly from Proposition |2J the minimum determinant of the infinite code is equal to 

^ ,^ , _ 1 .. ^ ^ _ 45 _ 1 _ 1 
5™„(Coo) - Y^ • NK,Q{a) _ _ _ ^^ _ ^. 

VI. 3x3 PERFECT STBC CONSTRUCTION 

In this case we use HEX symbols. Thus, the base field is F = Q(j). Let 6 = Cr + Cf ^ = 2cos(^) and 
K be Q(j, 9), the compositum of F and Q(6l). Since (/?(7) = 6, [Q(6') : Q] = 3, and thus [Q(i, 6) : F] = 3. 
The discriminant of Q(0) is di^tQ\ = 49 the minimal polynomial pe{X) = X^ + X"^ — 2X — 1. The 



extension K/F is cyclic with generator t : Ct + Qj "^ ^^ Qj + Cj 



2 

The corresponding cyclic algebra of degree 3 is ^ = {K/F,a,^), that is 



with e € A such that e^ = 7 G F* and le = ea{l) for all / G i^. In order to obtain a perfect code, we 
choose 7 = j. 
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A. The Z[j]-lattice Z[jf 

In this case, we look for a Z[j]-lattice which is a rotated Z[j]^{= A^) lattice. The relative discriminant 
of K is dx/F = d.Q/g\ = 49 = 7^, while its absolute discriminant is dx = —3^ • 7^. A necessary 
condition to obtain a rotated Z[j]^ lattice is the existence of an ideal T C Ok with norm 7. In fact, 
the lattice A{Ok) has fundamental volume equal to 2^'^y^pxl = 7^ f^J and the sublattice A(T) 

has fundamental volume equals to 2^^y^JdK\N{I) = 7'^ ( ^ j , where the norm of the ideal N{Z) is 
equal to the sublattice index. This suggests that the algebraic lattice A(T) could be a homothetic (scaled 
rotated) version of A2, namely, (7^42) . 

An ideal I of norm 7 can be found from the following ideal factorizations 

{7)Ok = I^It". 

Let us consider T = Tj. It is a principal ideal X = (a) generated by a = (1 + j) + 6. K Z[j]-basis of 
{a)OK is given by {ae^]l^^^ = {(1 + j) +0,{l+ j)9 + 9'^,! + 26 + jO'^}. Using the change of basis 
given by the following matrix 



1 











-1 


1 


2 


1 






TrK/Qij){l + 9 


+ 0^ 


) = 7 


if 


k = 


-- 1 


Trx/(Q(j)(2- 


-9)- 


= 7 


if 


k = 


-2 


. Tr^/Q(j)(4- 


0') 


= 7 


if 


k = 


-3 



one gets a reduced Z[j]-basis 

{j^fclLi = {(1 + j) + 0, (-1 - 2j) + j9^ (-1 - 2j) + (1 + j)9 + (1 + j)9''}. 
Then by straightforward computation we find 

-j^^K/Qij) i^m) = 5ki k, I = 1,2,3 

using TrQ(5i)/Q(l) = 3, TrQ(^QyQ{9) = -1, TrQ(^gyQ{6 ) = 5. 
We compute, for example, the diagonal coefficients 

Trii-/Q(j)(j^fej7fc) 

The generator matrix of the lattice in its numerical form is thus given by 
1 

V7' 

0.66030 + 0.32733i 0.02077 + 0.32733i -0.49209 + 0.32733i 

-0.29386 - 0.14567i -0.03743 - 0.58982i -0.61362 + 0.40817i 

0.52952 + 0.26250i -0.04667 - 0.73550i 0.27309 - 0.18165i 

B. The norm condition 

We show that the rank criterion is fuUfilled by this new code. The following proposition guarantees 
that A = {Q{j,9)/F,a,j) is a division algebra. 

Proposition 8: The units j and j"^ are not norms in Q{j,9)/F. 

Proof: See Appendix Hill for the proof, which uses Class Field Theory. ■ 
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C. The minimum determinant 

As the ideal T is principal, we can use dlOt or Proposition |2l to get 



1 



7 



1 



1 
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^mm(Coo) = -^Nx/q{ol) 



VII. 6x6 PERFECT STBC CONSTRUCTION 

As in the 3 antennas case, we transmit HEX symbols. Thus, the base field is F = Q(j). Let = C28 + 

C28^ = 2cos(yj) and K be Q(j, 9), the compositum of F and Q((9). Since c^(28) = 12, [Q(6») : Q] = 6, 
and thus [Q(j, 0) : F] = 6. The extension K/F is cyclic with generator cr : C28 + C28^ ""^ Cis + C28^- 
The corresponding cyclic algebra of degree 6 is ^ = (ET/F, cr, 7), that is 

A = l-K®e-K®e^-K®e^-K®e'^-K®e'-K 

with e ^ A such that e^ = 7 G F* and /e = ea{l) for all / G i^. In order to obtain a perfect code, we 
choose 7 = — J. 



A. r/ie Z\j]-lattice Z[jf 

First note that the discriminant of K is dx 



2^2 . 36 . 710 following the approach given in Section 



IIII-BI we need to construct a 



'[J\ 



lattice. 



A necessary condition to obtain a rotated version of Z[j]^ is that there exists an ideal T C Ok with 



v^V 



and 



norm 7. In fact, the lattice A{Ok) has fundamental volume equal to 2 ^\/\dK\ = ^^ ' 2^ ■ 1 -5- 1 

the sublattice A(T) has fundamental volume equal to 2^^^/\dK\N{I) = 7'' • 2^ • (^j , where the norm 
of the ideal N{I) is equal to the sublattice index. This suggests that the algebraic lattice A(X) could be 
a homothetic version of A^, namely, (\/l4A2) , but this needs to be checked explicitly. 
An ideal I of norm 7 can be found from the following ideal factorizations 

-6 



(7)0 



K 



I J If 



Let us consider 1 = X-j. Unlike in the preceeding constructions, the ideal X is not principal. This makes 
harder the explicit computation of an ideal basis, and in particular of the ideal basis (if any) for which 
the Gram matrix becomes the identity. 

We thus adopt the following alternative approach. We compute numerically a basis of X, from which 
we compute a Gram matrix of the lattice. We then perform a basis reduction on the Gram matrix, using 
an LLL reduction algorithm (see Appendix IVII for more details). This gives both the Gram matrix in the 
reduced basis and the matrix of change of basis. We get the following change of basis 
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and the lattice generator matrix in numerical form 

/ 1.9498 1.3019 -0.8660i -0.0549- 0.8660i -1.7469 - 0.8660i 1.5636 0.8677 \ 

0.8677 -1.7469 -0.8660i 1.3019-0.8660? -0.0549 - 0.8660i -1.9498 1.5636 

1.5636 -0.0549 -0.8660i -1.7469-0.8660? 1.3019 - 0.8660i -0.8677 -1.9498 

-1.9498 1.3019 -0.8660i -0.0549- 0.8660i -1.7469 - 0.8660i -1.5636 -0.8677 

-0.8677 -1.7469 -0.8660i 1.3019 - 0.8660i -0.0549 - 0.8660i 1.9498 -1.5636 

V -1.5636 -0.0549 -0.8660i -1.7469- 0.8660i 1.3019 - 0.8660i 0.8677 1.9498 / 

This matrix satisfies RR^ is the identity matrix, so that we indeed get a rotated version of the Ag 
lattice. 

B. The norm condition 

Since 7 = —j, we have to check that —j,j'^, —j^ = —l,j* = j and —j^ = — j^ are not norms in K. 
Lemma 3: We have the following field extensions: 

Q(j) c Q(j, C7 + C7"') c Qij, C28 + Css') 

Proof: The proof is similar to that of Lemma |2l One has to show that Q(i, C7 + Cf^) is the subfield 
fixed by (a^), the subgroup of order 2 of Gal(Q(j, C28 + C28')/F) = {a). ■ 

Proposition 9: The algebra A = {K/F, a, —j) is a division algebra. 

Proof: We prove, by contradiction, that ±j and itj^ are not norms in K*. Suppose that either ±j 
or ibj^ are a norm in K*, i.e., there exists x G K* such that Nj^ip{x) = itj (resp. ±j^). By Lemma |3] 
and transitivity of the norm, we have 

^K/F{x)=NQ^^-^^^^^-iyp{Nj^/Q^^-^^^^^-i^{x)) = ±j (resp. ±f). (20) 

Thus j and j"^ have to be a norm in Q(j, (,7 + C?^^), which is not the case, by Propositions 1 1 01 and 1 1 1 1 in 
Appendix ITO] 

For the cases of —j and — j^, since [Q(j, C7 + Cf ^) : ^] = 3, d20l yields 

^Qa,c.+c.-^)/F(-^i^/Q(i,c.+c.-^)(^)) = i (""^^P- ^■')' 

which gives the same contradiction. 

The proof that - 1 is not a norm can be found in Appendix |V] and uses Class Field Theory. ■ 

C. The minimum determinant 

Since the ideal X is not principal, we use the bounds of Corollary |3l 

!,,.„. 1 1 „ .. . 1 .., . 72 



TTfi • ^kU^) = T^fT^ = 1 < <5™„(C^) < —. minN{x) 



^46 ^-A/^v-; 26 . 75 dQ(,) ^ "™"^-°°^ - 146 Tex ^■-' 26 • 76 

yielding 

1 1 



26 . 75 - "mmv-oo; ^ 26 • 74 ■ 
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VIII. Existence of perfect codes 

Since we have given constructions only for dimensions 2, 3, 4 and 6, it is interesting to discuss the 
existence of perfect codes. Perfect space-time block codes must satisfy a large number of constraints. Let 
us derive here the consequences of these constraints in the choice of the corresponding cyclic algebra. 

First note that in order to have non vanishing determinants when the spectral efficiency increases, 
determinants of the infinite code C^o must take values in a discrete subset of C. We have shown in 
Section UlI-CI that the determinants of Cj are in Op, when Z C Ok and 7 E Op- But Op is discrete 
in C if and only if F is a quadratic imaginary field, namely F = Q{\/—d), with d a positive square 
free integer. Indeed, we have that \a + 6\/^p G Z if a, 6 G Z. The positive minimum of an integer is 
thus 1. This is not true anymore if we consider already \a + 6\/dp, which belongs to Z[\Q]. We cannot 
obtain a minimum without any constraint on a, 6 G Z. The same phenomenon appears even more clearly 
in higher dimension. 

The average energy per antenna constraint requires I7I = 1. Furthermore, the proof of the non-vanishing 
determinant relied on 7 being in Op- There are two ways of getting a tradeoff between these two 
conditions. Our approach consists in choosing 7 to be a root of unity. Since the base field has to be 
quadratic, this gives as choice Q{i), which contains the 4th root of unity i, and Q(j), which contains 
the 3rd root of unity j and the 6th root of unity —j. The following lemma confirms these are the only 
possibilites: 

Lemma 4: [25, p.76] Let d be a positive square free integer. The only units of F = Q{V—d) are ±1 
unless F = Q(i) or F = Q(j). 
As a consequence, the perfect codes proposed are available only in dimension 2, 3, 4 and 6. 

Elia et al. recently considered the option of droping one of the two conditions. In [13], they drop the 
constraint I7I = 1, at the price of loosing the average energy advantage. They also consider an element 
7 of norm 1, but not in Op. Since 7 = 71/72 G F, the minimum determinant of the resulting code can 
be written as -. — m^rrjj det(X), where X is a codeword with coefficients in Ok- Thus the non-vanishing 
determinant property holds, but there is a loss in the coding gain proportional to |72n'^~^^ These codes 
are not restricted to the dimensions 2, 3, 4 and 6. 

IX. Simulation results 

We have simulated the complete MIMO transmission scheme using perfect Space-Time codes, and the 
previously best known codes. Transmitted symbols belong to q-QAM (2 and 4 antennas) or g-HEX (3 
antennas) constellations, q = 4,8, 16, 64. We used the modified version of the Sphere-decoder presented 
in [29]. 

QAM constellations have minimum Euclidean distance 2. The respective average energy per symbol for 
the 4, 8, 16 and 64-QAM constellations are 2, 6, 10 and 42. The q-HEX constellations are finite subsets of 
the hexagonal lattice A2. In fact the hexagonal lattice is the densest lattice in dimension 2; constellations 
using points from the hexagonal lattice ought to be the most efficient [14]. Since ^42 is not a binary lattice, 
bit labeling and constellation shaping must be performed ad hoc. The best finite hexagonal packings for 
the desired sizes are presented in Figures |2l |3l and |3 
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The respective average energy per symbol for the 4, 8, and 16 HEX constellations with minimum 
Euclidean distance 2 are 2, 4.5 and 8.75. We should note the energy saving compared to QAMs of the 
same size. The HEX constellations are carved from shifted versions of the lattice 2A2. For 4, 8-HEX 
constellations the shift is in (1,0), while for 16-HEX constellation the shift is in (1/2,0). 

In Figure |51 we have plotted the codeword error rates for the Golden code (GC), some other 2x2 
Perfect codes (PC) and the best previously known 2x2 STBCs [9] (BPC), as a function of E^j/Nq, using 
4, 16,64-QAM constellations. In [9], the values of 7 giving the best codes were obtained by numerical 
optimizations and depend on the spectral efficiency. As we concluded in [4], [5], the Golden code has 
the best performance. We see in Fig. |5l that perfect codes with p = \/l3 and p = \/37 have performance 
close to that of the BPCs. However the code with p = \/T7 which is not a perfect code (the cyclic 
algebra is not a division algebra) has the worst performances, and we can even observe a change in the 
slope of the curve for high SNR, due to the reduced diversity order of this code (2 instead of 4). In fact, 
as shown in Example |3 there exists an x ^ K = Q(i,\/l7) such that Nx/Q(i){x) = i. The appearance 
of such an x is rare, which explains why this code works well at low and medium signal to noise ratio 
and the change of slope appears at very low error rates. 

In Fig. |6l and Q we have plotted respectively the codeword error rates of the 3x3 and the 4 x 4 PC 
and the best previously known codes [11], [16] as a function of Eb/No. In Fig. |6l we see that for the 
4-HEX constellation, the BPC performs a little better than the PC. However, when the constellation is 
8-HEX or larger, PCs have better performance, due to the constant minimum determinant. 

In Fig. Q we note that the 4x4 PC improves over BPC codes when we use the 64-QAM. 

X. Conclusion 

In this paper we presented new algebraic constructions of full-rate, fully diverse 2x2, 3x3, 4x4 
and 6x6 space-time codes, having a constant minimum determinant as the spectral efficiency increases. 
The name perfect STBC, used for these codes, was suggested by the fact that they satisfy a large number 
of design criteria and only appear in a few special cases as the classical perfect error correcting codes, 
achieving the Hamming sphere packing bound. 



Appendix I 
Number Fields: Basic Definitions 

The codebooks we build are based on cyclic algebras built over number fields, we thus need some 
background on number fields. This appendix aims at giving intuition to the reader who does not know 
the topic. It focuses on examples, and may skip some technical points in order to be more accessible. 

Number fields can first be thought of as finite vector spaces over a base field. For example, Q{i) = 
{a + bi, a,b ^ Q} is a vector space of dimension 2 over Q, whose basis is given by {1, i}. In our case, 
we will consider two number fields, denoted by K and F, and K will be a vector space of dimension 
n over F. We say that i^ is a field extension of F, which we denote by K/F. The dimension of K 
over F as a vector space is called the degree, and is denoted by [K : F]. Another way of thinking of 
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a number field is to add a root of a polynomial, with coefficients in F, to a field, and to add also all 
its powers and multiples, so that the resulting set is indeed a field. For example, Q(i) is built adding 
the roots of the polynomial X"^ + 1 to Q. The field extension K/F can similarly be seen as adding the 
element 9, root of a polynomial p{X), to K. We may write K = F{6). Since a polynomial has n roots, 
one may wonder if taking one root or another may change the number field. If all the roots are indeed in 
the number field, it does not change, and the number field is called a Galois extension. Not all number 
fields are Galois extensions. 

For our purpose, we are interested in a field extension K/F such that all roots 9i,...,9n of p{X) 
are not only in K, but furthermore are related to each other as follows: there exists a map a such 
that a^{9i) = 9j, k,j = l,...,n. In such case, K/F is called a cyclic Galois extension, and {cr^}, 
k = I, . . . ,n, is called a (cyclic) Galois group (it can be shown that it has indeed a group structure). For 
example, Q(i) is a cyclic Galois extension of degree 2, since there exists a : i \-^ —i. 

There are two important objects that can be defined thanks to {<7^}, k = 1, . . . ,n. We define the trace 
and the norm of an element x G K resp. as follows: 

n— 1 n— 1 

k=0 fe=0 

We may call a relative trace/norm if the base field is not Q, by opposition to an absolute trace/norm 
when F = Q. 

Let now L be a number field of degree n over Q. Consider the set of elements x of L that satisfy the 
following property: there exists a monic polynomial / with coefficients in Z such that f{x) = 0. This 
set is called the ring of integers of L, and is denoted by Ol- It can be shown that this is indeed a ring, 
but what is more interesting is that this set posseses a Z-basis. We will use this fact extensively in the 
paper. Let {uji, . . . ,u;„} be a Z-basis, i.e., all elements can be written as integer linear combinations of 
basis elements. Then det{TY ];^/Q{LOiUJj)f j^^) is an invariant of L called the discriminant of L. Similarly 
as for the trace and norm, we call the discriminant absolute to emphasize that the base field is Q, and 
relative otherwise. 

Appendix II 
The Hasse Norm Symbol 

In this Appendix, we introduce the Hasse Norm Symbol. It is a tool derived from Class Field Theory, 
that allows to compute whether a given element is a norm. Our exposition is based on [18]. In the 
following, we consider extensions of number fields K/F that we assume abelian. 

Denote by K^ the completion of K with respect to the valuation v. We denote the embedding of K 
into K^ by i^. 

Definition 5: [18, p. 105] Let K/F be an abelian extension of number fields with Galois group 
Gal{K/F). The map 



^) : K* -^ Ga.\{K/F) 



i^(x),K/F 

V 
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is called the Hasse norm symbol. 

The main property of this symbol is that it gives a way to compute whether an element is a local norm 

[18, p. 106,107]. 

Theorem 2: We have ( ^' J ) = 1 if and only if x is a local norm at v for KjF . 
In order to compute the Hasse norm symbol, we need to know some of its properties. Let us begin with 
a property of Unearity. 

Theorem 3: We have 

' xy,K/F \ ^ / x,K/F \ f y,K/F \ 

We then know how the symbol behaves at unramified places [18, p. 106]. 
Theorem 4: If v is unramified in K/F, then we have, for all x £ F*: 

x,K/F \ _ /K/FV^^) 

where ( — ^ 1 denotes the Frobenius of v for K/F (see Remark|3]below), and v{x) denotes the valuation 
of X. 

Remark 3: For our purpose, it is enough to know that the Frobenius ( — ^ J is an element of the 
Galois group Ga[{K/F). We do not need to know it explicitly. For a precise definition, we let the reader 
refer to [18, p. 107]. 

Corollary 6: At an unramified place, a unit is always a norm. 

Proof: It is straightforward since the valuation of a unit is 0. ■ 

A remarkable property of the Hasse norm symbol is the product formula [18, p. 113]. 

Theorem 5: Let K/F be a finite extension. For any x G F* we have: 

n(^)-. 

where the product is defined over all places v. 

Remark 4: By Corollary |6l we know that a unit is always a norm locally if the place is unramified. 
Since we are interested in showing that a unit 7 is not a norm, we will look for a contradiction at a 
ramified place. 

Before giving the proofs in themselves, we explain briefly their general scheme. The idea is to start 
from the product formula, and to simplify all the terms except two in the product over all primes, so that 
we get a product of two terms equal to 1 : 

Hopefully, one of the two terms left will involve 7, the other will be shown to be different from 1, so 
that since the product is 1 , we will deduce that the term involving 7 is different from 1 , thus 7 is not a 
norm. In order to make it easier to simplify the product formula, we introduce an element y ^ K such 
that y7 is a unit locally at ramified primes, and we compute the product formula 

n (^^) - - 
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Appendix III 



j AND j2 ARE NOT A NORM IN Q{j, 2 COS '' — 



7 

In this section, we prove that j and j^ are a not a norm in Q{j, 2 cos (^) )/Q(j). We show that j and 
j"^ are not a norm locally by computing their Hasse norm symbol. The proof is detailed for j. 
Proposition 10: The unit j is not a norm in K/F = Q(j, 2cos (^))/Q(i)- 
Proof: We consider the field extension K/F. We have 

7Z[j] = (i-2)(j + 3)=P7q7. 
We show that j is not a norm locally in pj, thus j is not a norm in K. 

We look for a number y in Z[j] satisfying 

y = 1 (mod j - 2) (21) 

jy = 1 (modj + 3). (22) 

By applying the Chinese Remainder Theorem over Z[j], we find y = 7 — 3j with (y)Z[j] = pyg. Let 



- j denote the Hasse norm symbol. By the product formula 

n(^)- n (^) n (^)-- 

^ i/ ramified ;/ unramified 

The product on the ramified primes yields ( -^^' ' j ( -'^' ^ j, since the ramification in K/Q{j) is 

in 7 only. Note that ( ^^\ ) = ( ^' ^ ) ( ^' ^ j by linearity. We now look at the product on the 
unramified primes. Since y G pyg, its valuation is zero for v ^ pyg. The valuation of a unit is zero for 
all places, so that we get 

n (^)^ n (^)(^)K^)- 

u unramified u unramified 

Thus equation \23\ simplifies to 

'3,K/F\ fy,K/F\ fjy,K/F\ fy,K/F\ 



P7 J \ P7 J \ (]7 J \ P79 J 

The second and third terms are 1 by choice of y (see equations d2lT) and d22b '). so that finally we have 

1,K/F\ f y,K/F \ ^ ^ 

p7 y V p79 y 

Since P79 is inert, the second term is different from 1, so that ( ■'' ' j / 1. In words, j is not a norm 
in p7 which concludes the proof. ■ 

Proposition 11: The unit f is not a norm in K/F = Q(j, 2cos (^))/Q(j)- 

Proof: The proof that j"^ is not a norm is similar to the above one. We keep the notation of the 
above proof. We show that j^ is not a norm locally in p7, thus j^ is not a norm in K. 
Let y = 5j — 9. We have that 

y = 1 (mod j - 2) (24) 

fy = 1 (mod3 + j) (25) 
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and {y)'L[j] = pisi. Repeating the same computations as in the above proof, we get 

( j\K/F \ f y,K/F \ ^ ^ 

V P7 y V Pi5i ) 
where pisi is inert. This imphes that j"^ is not a norm. ■ 



Appendix IV 

15 



i^ = -1 IS NOT A NORM IN Q(i, 2 COs(^ 



We prove here that i^ = — 1 is not a norm in Q(i, 2cos(||-))/Q(i). The general scheme of the proof 
is the same as in Appendix ITO] though we have to be a bit more careful here, since the ramification in 
Q(i,2cos(f|))/Q(i) appears in two primes, unlike in Q(j,2cos (^))/Q(j). 

Proposition 12: The unit -1 is not a norm in K/F = Q(i, 2cos(f|))/Q(i). 
Proof: We consider the field extension K/F. We have 

5Z[i] = {i + 2){i - 2) = pgqs and 3Z[i] = 3 = ps. 

We show that i is not a norm locally in ps, thus i is not a norm in K. We look for a number y in Z[i] 
satisfying 

y = 1 (mod i + 2) (26) 

-y = 1 (mod i - 2) (27) 

-y = 1 (mod 3) (28) 



By applying the Chinese Remainder Theorem over Z[i], we find y = 12i — 25 with {y)Z[j] = pyeg. Let 
( ^' J ] denote the Hasse norm symbol. By the product formula 

n (=^) - '■ 

The product on the ramified primes yields ( ~^' ' j ( ~^' ' ) ( ~^' ' ) , since the ramification in 
KjF is only in 5 and 3. Since y G p769, its valuation is zero for v ^ p769- The valuation of a unit is 
zero for all places, so that we get for the product on the unramified primes 

n {^^y n (-^)(^)^(^)- 

V unramified v unramified 

Thus equation (I29t simplifies to 

( -y,KlF \ ( y,K/F \ / -1,K/F \ / -y,K/F \ ( y,K/F \ ^ ^ 

V p3 y V p5 y V p5 y V is y V p769 y 

The first, second and fourth terms are 1 by choice of y (see equations d26b . d27t and ( I28t '). so that finally 
we have 



f -l,K/F \ / y,K/F \_^^ 

V p5 y V p769 y 



Since pygg does not split completely, the second term is different from 1, so that f — in^) v^ -*-' which 
concludes the proof. 
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Appendix V 

(-j)3 = -1 IS NOT A NORM IN Q(j, C28 + (2^8^ 

We prove here that (— j)^ = —1 is not a norm in Q(j, ^28 + C28^)/Q(i)- The proof is similar to that 
of Appendix 11 VI 

Proposition 13: The unit -1 is not a norm in K/F = Q(C28 + C28 ; i)/Q(i)- 
Proof: We consider the field extension K/F. We have 

7Z[j] = (i - 2)(j + 3) = P7q7 and 2Z[j] = 2 = pa- 

We show that —1 is not a norm locally in p7, thus —1 is not a norm in K. 
We look for a number y in Z[j] satisfying 

y = 1 (mod j - 2) (30) 

-y = 1 (mod3+i) (31) 

-y = 1 (mod 2) (32) 

By applying the Chinese Remainder Theorem over Z[j], we find y = 3 — 8j with (y)Z[j] = P97. Let 
( ^' J j denote the Hasse norm symbol. By the product formula 

n(^^)-. 

The product on the ramified primes yields ( ~^' ^ j ( ~^' ' j ( ~^' '^ j , since the ramification in 
K/F is in 7 and 2 only. Since y € P97, its valuation is zero for v ^ P97. The valuation of a unit is zero 
for all places, so that we get for the product on the unramified primes 

n {^^)- n (-^)(^)^(^)- 

V unramified v unramified 

Thus equation (I33t simplifies to 

-y,K/F \ ( y,K/F \ / -1,K/F \ / -y,K/F \ ( y,K/F \ ^ ^ 
p2 J \ P7 J \ P7 J \ q? J \ P97 J 

The first, second and fourth terms are 1 by choice of y (see equations d30b . d^TT i and (I32t ). so that finally 
we have 

-1,K/F\ fy,K/F\ 



(^) 



P97 / 



1. 



Since P97 does not split completely, the second term is different from 1, so that ( — '—^ j 7^ 1, which 
concludes the proof. 
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Appendix VI 
The LLL reduction algorithm over Z[j] 

The standard LLL reduction algorithm [19] over Z, can be easily modified to work over Z[j] [23]. The 
two main points to be careful about are 

• the Euclidean division: the quotient of the Euclidean division over Z[j] is defined as follows: let 
X = xi + jx2 and y = yi + jy2, xi,X2,yi,y2 G ^- The division of x by y yields | = -^i + JZ2, 
with zi,Z2 G Q. Then we have that x = yq + r, where q = [zi] + j[z2\- 

• the conjugation: the usual complex conjugation is replaced by the r-conjugation, that sends j onto 
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Fig. 4. 16-HEX Constellation. 
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Fig. 5. Golden code (GC) (perfect code with p = 5) and other perfect codes (PC) compared to the best peviously known 
2x2 codes (BPC). 
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Fig. 6. Perfect 3x3 codes with HEX symbols (PC) compared to thie best peviously known codes (BPC). 



DRAFT 



February 1, 2008 



OGGIER ET AL. PERFECT STBCS 



37 




-A PC, 4-QAM 
-V PC, 16-QAM 
-X PC, 64-QAM 
G — O BPC, 4-QAM 
D — a BPC, 16-QAM 
O — O BPC, 64-QAM 



15 

Eb/NO (dB) 



Fig. 7. Perfect 4x4 codes with QAM symbols (PC) compared to the best peviously known codes (BPC). 
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